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Sphaleron transition rate in the classical 1 + 1 dimensional 
abelian Higgs model at finite temperature 
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We compute the sphaleron transition rate in the 1 + 1 dimensional abelian Higgs model at finite temperature, 
by real time simulation using the classical canonical ensemble. 



1. INTRODUCTION 

A nonperturbative calculation of the sphaleron 
transition rate in the Standard Model at finite 
temperature is difficult because it involves real 
time correlations. To deal with this the clas- 
sical approximation was proposed and tested in 
the 1+1 dimensional abelian Higgs model [1], 
for which a semiclassical calculation of the rate 
was available [2]. The test involved the mi- 
crocanonical ensemble and sphaleron transitions 
were counted 'by hand'. A nonperturbative defi- 
nition of the rate 

A{t) = {[Cit) - C7(0)]2) = R, large t, (1) 

where C'(t) the Chern-Simons number, was used 
in a real time Langevin simulation [3]. This re- 
produced the temperature dependence and may 
be useful for a determination of the absolute rate 
[4]. Interesting results have also been obtained 
recently with a realistic heat bath method [5]. 

Here we shall present results using the canoni- 
cal ensemble, 

A(t) = J DpDqeM-H{p,q)/T] 

[Cipit),qit))-Cip,q)r. (2) 

The canonical variables p(t) and q(t) are evalu- 
ated with Hamilton's equations from the initial 
conditions p(0) = p, q(0) = q, with p,q drawn 
from the ensemble. 

The effective action of the abelian Higgs model 
at finite temperature will have approximately the 
classical form Seff ~ -/ dt dx [Ff.^Pt"' /Ag^ + 
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(d^<i> - iA^<i>)*(d''<i> - iA"<j>) + m^\<j>\^ + \\<j>\% at 
weak coupling A/m^ ^ 1, g^/\ = 0(1), with 
temperature dependent parameters m^. A, . 
This (quantum) effective action is used for the 
classical approximation. 

In the Higgs phase < 0, = — 2m^, m\ = 
((/^/2A)m^, and the expectation value {(j>) = v/2, 
serves an expansion parameter = 2A/m^. 

The sphaleron energy is given by Eg = 2v^m^/3. 

At high temperatures T/mA,,p ^ 1, we may 
find classical behavior. Keeping the ratio Eg/T 
fixed this means ^ 1, i.e. the semiclassical 
region. 

To deal with the Rayleigh-Jeans divergence the 
classical field theory is put on a lattice, with lat- 
tice distance a of order of the inverse temperature, 
a = c/T, with c a numerical constant. It is con- 
venient to make a scale transformation to dimen- 
sionless variables [3], a ^ a 
(f) ^ vcf), v\/^Afi, which gives an overall 

factor in the action. The scaled lattice action 
takes the form 

. N-l 

S = dtaY,{[{Atn+i-At^)/a-dtA,^f/2^ 

■' n = 

+ |(exp(-iaAj;„)(?i„+i - (?i„)/ap 
-|(5t-iAt„),/-„P + (|,/,„|2-l/2)2},(3) 

where ^ = (/^/A and we assume periodic bound- 
ary conditions with spatial size L = Na. The 
hamiltonian and sphaleron energy scale as H ^ 
v^\/\H , Eg v^\/\Es. It is convenient to scale 
also the temperature T v^\/XT, such that 
the Boltzmann factor keeps its usual form. In 
scaled variables {<j>) = 1/^2, Eg = 2^2/3 ^ 0.94, 



2 



= V2, niA = v^. 
A crucial observation is that the scaled lattice 
distance a goes to zero in the semiclassical limit: 
under the scaling we have c = aT c = aTv^ . 
Hence a oc at fixed scaled temperature (fixed 
Eg/T). In the classical approximation drops 
out, but its role is taken over by the scaled lattice 
distance. 

Following ref. [3] we use the Coulomb gauge 
dxAx = and solve for the Gauss constraint. 
Then the Chern-Simons number 



C 



dx Ax/2tt = LAx/2tt 



(4) 



can be taken as one of the canonical variables. 
The hamiltonian follows as usual. 

2. NUMERICAL STRATEGY 

We use a second order Langevin procedure to 
generate the canonical ensemble of p's and g's. 
Following ref. [3] we use polar coordinates (f) = 
pexTp(i6) and apply the random Langevin forces 
only to the gauge invariant variables p, pp, and 
not to the gauge variant 0, pg. The random forces 
are also not applied to C and pc , which makes it 
possible to monitor thermalization of pc- The 
system is kicked at times t = kh (integer k), and 
for H = X]n^'pn/4fl + • • • the algorithm can be 
presented as 



Pn(t + h) 



Pn(t) 



, , OH 

d^l, 

t opp^ 



PpJt) + ^AafhT^n{t) 



(5) 



(6) 



where / is the friction coefficient, T is the speci- 
fied ('input') temperature and the rj are indepen- 
dent gaussian random numbers with unit variance 
and zero mean. The rj are absent from the other 
dynamical equations. The singularity of polar co- 
ordinates at /9 = led us to use a variable stepsize 
leapfrog algorithm for the integration of Hamil- 
ton's equations over the interval (t,t + h). 

The finite Langevin stepsize h introduces an er- 
ror which we monitored by requiring the 'output 



temperatures' Tc and Tp, 
,CLpl 



47r2 



Tc 



'2a' 



P' 



(7) 



to be such that their inverse l3c^p = l/Tc^p dif- 
fered less than 0.1 in absolute value from /3 = 
l/T. This criterion is based on the interpretation 
that Langevin errors lead to an effective temper- 
ature, the output temperature. Then the expec- 
tation r cx exp(— /JE's), suggests that an absolute 
error in /3 of 0.1, causes a relative error in F of 
^ 10%. 

Having produced an independent configuration 
of p's and g's we took this as the initial condi- 
tion for the real time integration of Hamilton's 
equations. For this we could use the original 
cartesian coordinates, as the condition of charge 
zero is easily satisfied by projecting regularly onto 
zero charge (this only involves changes of machine 
precision order). After real time integration the 
Langevin process was started again and the pro- 
cess was repeated untill sufficient statistics was 
obtained. 

Our configurations were actually more micro- 
canonical than canonical, because the Langevin 
processes were stopped (for historical reasons) at 
times that the total energy had its mean value. 
As expected we found in a few checks that the 
true canonical ensemble gives the same results 
within errors. The friction parameter was taken 
/ = 1, and we checked that results do not depend 
on /, as expected. 
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Figure 1. Difl^usion A{t) for /3 = 11, ^ = 4, L = 
16, a = 0.32. 

For low temperatures A tends to be small and 



3 



large real times t were needed. We required A > 
10. An accurate fifth order Adams-Bashforth- 
Moulton predictor-corrector multi-step integra- 
tion algorithm was used to keep the numerical 
drift in the total energy sufficiently small. Oth- 
erwise the diffusion A will not be linear in t (cf. 
fig. 1). 

3. RESULTS 
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Figure 2. Winding number (line) and Chern- 
Simons number (dots) for /3 = 8, ^ = 4, _L = 16, 
and a = 0.32, during Hamilton evolution. 

Fig. 2 shows an example of the real time be- 
havior of C at low temperature. A comparison 
is made with the winding number of the scalar 
field w = - J2„ dO„, with d6„ = 6„+i - 0„ (mod) 
27r G (— 7r,7r]. On the lattice the winding and 
Chern-Simons numbers are defined modulo N . 
Fig. 1 shows an example of A(t) at low temper- 
ature. In fig. 3 a comparison is made with the 
analytic semiclassical result 



(3) 



1/2 



mlL 



(27r)3/: 



ipE.f' eM-PE,) (8) 
r(a-Fs-F l)r(a-s)l^''^ 



T{a+ l)r(a) 



where a = ^ and s = (-1 ^1 + 8^/2. The 
value of ^ is again ^ = 4. The upper data are 
for a = 0.32, the lower data for a = 0.16 (at 
l3 = 10, 11 for a = 0.32 only). The system size is 
L = 16 (# = 50 and N = 100), and we obtained 



the same results within errors for L = 10.28, for 
a number of /3 values. However, for sizes as low 
as _L = 8 we see clear deviations. 
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3. Results for ln(r/m^i.). 



The solid line 



represents the analytic formula (8). 

The errors in the rate F are statistical and de- 
termined with the jackknife method from at least 
900 configurations. For /3 we used the output 
temperature (7), with errors obtained by the bin- 
ning method. The input /3 is the nearest integer. 

We conclude from fig. 3 that the classical sim- 
ulation is able to reproduce the semi-classical for- 
mula for /3^7. The data at smaller (] depend on 
the lattice distance a and appear therefore to be 
ambiguous. Perhaps smaller a are needed in this 
region of /3 values; or the classical or the semi- 
classical approximation or both break down. 
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